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ABSTRACT 

We present constraints on the mean dark energy density, fix and dark energy equa- 
tion of state parameter, u>x, based on Chandra measurements of the X-ray gas mass 
fraction in 26 X-ray luminous, dynamically relaxed galaxy clusters spanning the red- 
shift range 0.07 < z < 0.9. Under the assumption that the X-ray gas mass fraction 
measured within r25oo is constant with redshift and using only weak priors on the 
Hubble constant and mean baryon density of the Universe, we obtain a clear de- 
tection of the effects of dark energy on the distances to the clusters, confirming (at 
comparable significance) previous results from Type la supernovae studies. For a stan- 
dard ACDM cosmology with the curvature f^K included as a free parameter, we find 



Q\ = 0.94 



+ 0.21 

-0.23 



per cent confidence limits). We also examine extended XCDM 



dark energy models. Combining the Chandra data with independent constraints from 
cosmic microwave background experiments, we find fix — 0.75 ±0.04, fl m = 0.261q'q4 
and wx = — 1.26±0.24. Imposing the prior constraint wx > — 1, the same data require 
wx < —0.7 at 95 per cent confidence. Similar results on the mean matter density and 



dark energy equation of state parameter, Q m = 0.24 ± 0.04 and wx = —1.20fo^g, are 
obtained by replacing the CMB data with standard priors on the Hubble constant and 
mean baryon density and assuming a flat geometry. 

Key words: X-rays: galaxies: clusters - cosmological parameters - dark matter - 
cosmic microwave - gravitational lensing 



1 INTRODUCTION 

The matter content of the largest clusters of galaxies is 
thought to provide an almost fair sample of the matter con- 
tent of the Universe (e.g. White et al. 1993; Eke et al. 1998). 
The observed ratio of baryonic-to-total mass in such clusters 
should therefore closely match the ratio of the cosmological 
parameters fib /fim, where fib and fi m are the mean baryon 
and total mass densities of the Universe in units of the criti- 
cal density. The combination of robust measurements of the 
baryonic mass fraction in the largest galaxy clusters with 
accurate determinations of fib/i 2 from cosmic nucleosynthe- 
sis calculations (constrained by the observed abundances of 
light elements at high redshifts) and/or cosmic microwave 
background (CMB) experiments, with a reliable measure- 
ment of the Hubble constant, Hq, can therefore be used to 
determine fi m . 

This method for measuring fi m , which is particularly 
attractive in terms of the simplicity of its underlying as- 



sumptions, was first highlighted by White & Frenk (1991) 
and subsequently employed by a number of groups (e.g. 
Fabian 1991, White et al. 1993, David, Jones & Forman 
1995; White & Fabian 1995; Evrard 1997; Mohr, Mathiesen 
& Evrard 1999; Ettori & Fabian 1999; Roussel, Sadat & 
Blanchard 2000; Allen, Schmidt & Fabian 2002a; Allen et 
al. 2003; Ettori, Tozzi & Rosati 2003; Sanderson & Ponman 
2003; Lin, Mohr & Stanford 2003). In general, these studies 
have found fi m < 1 at high significance, with recent work 



favouring fi n 



0.3 hi 



Sasaki (1996) and Pen (1997) were the first to describe 
how measurements of the apparent redshift dependence of 
the baryonic mass fraction could also, in principle, be used to 
constrain the geometry and, therefore, dark energy density 
of the Universe. The geometrical constraint arises from the 
dependence of the measured baryonic mass fraction values 
on the assumed angular diameter distances to the clusters. 
Although theory and cosmological simulations suggest that 
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the baryonic mass fraction in the largest clusters should be 
invariant with redshift {e.g. Eke et al. 1998), this will only 
appear to be the case when the reference cosmology used in 
making the baryonic mass fraction measurements matches 
the true, underlying cosmology. The first successful appli- 
cation of such a test was carried out by Allen, Schmidt & 
Fabian (2002a; hereafter ASF02; see also Allen et al. 2003a) 
using a small sample of X-ray luminous, dynamically re- 
laxed clusters with precise mass measurements, spanning 
the redshift range 0.1 < z < 0.5. These authors found 
O m = 0.30 ± 0.04 and ft A = 0.95t£^. A similar analysis 
was later carried out by Ettori, Tozzi & Rosati (2003) who 
obtained consistent results using a larger cluster sample that 
extended to higher redshift, but which included clusters with 
a wider range of dynamical states. 

The baryonic mass content of galaxy clusters is dom- 
inated by the X-ray emitting intracluster gas, the mass of 
which exceeds the mass of optically luminous material by a 
factor ~ 6 (e.g. White et al. 1993; Fukugita, Hogan & Pee- 
bles 1998; other mass components in clusters are expected 
to make only very small contributions to the total baryon 
budget). Since the emissivity of the X-ray emitting gas is 
proportional to the square of its density, the gas mass pro- 
file in a cluster can be precisely determined from X-ray data. 
Measuring the total mass profile is more difficult, however, 
and requires both a direct measurement of the X-ray gas 
temperature profile and the assumption of hydrostatic equi- 
librium in the gas. Measurements of the temperature profiles 
in intermediate-to-high redshift galaxy clusters have only be- 
come possible following the launch of the Chandra X-ray 
Observatory. The exquisite spatial resolution of Chandra 
makes measuring the temperature profiles of even distant 
clusters a relatively straightforward task, given sufficient ex- 
posure time. The use of the hydrostatic assumption in mak- 
ing the mass and baryonic mass fraction measurements is 
more problematic, however, and requires a restriction to dy- 
namically relaxed systems when carrying out a cosmological 
test similar to that described here. 

In this paper we present a significant extension to the 
ASF02 work. The cluster sample is significantly larger and 
includes 26 X-ray luminous, dynamically relaxed systems 
spanning the redshift range 0.07 < z < 0.9. As well as en- 
hancing the sample, we have also expanded the analysis. We 
now include a bias factor, b, in the X-ray analysis (see also 
Allen, Schmidt & Bridle 2003) , motivated by gas-dynamical 
simulations, that accounts for the (relatively small amount 
of) baryonic material expelled from such clusters as they 
form. We also examine XCDM dark energy models in which 
the equation of state parameter, t»x, is allowed to take any 
constant value. Finally, as well as results based on the Chan- 
dra data using simple priors on Qbh 2 , h and b, we also 
present results from the combination of Chandra and CMB 
data. This latter combination is shown to be particularly 
powerful in constraining the overall dark energy density. 

As in ASF02, we report measurements of the X-ray gas 
mass fraction for two reference cosmologies: an SCDM cos- 
mology with h = ffo/lOOkms -1 Mpc _1 = 0.5, n m = 1 and 
Qa = 0, and a ACDM cosmology with h = 0.7, fi m = 0.3 
and ft A = 0.7. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample selection 

Our sample consists of 26 galaxy clusters spanning the red- 
shift range 0.07 < z < 0.9, with X-ray temperatures kT > 
5keV and X-ray luminosities Lx, 0.1-2.4 > 10 45 ergs -1 . 
The clusters exhibit a high degree of dynamical relaxation in 
their Chandra images (sharp central X-ray surface bright- 
ness peaks, regular X-ray isophotes and minimal isophote 
centroid variations) and show no evidence for a signifi- 
cant loss of hydrostatic equilibrium in X-ray pressure maps 
and/or gravitational lensing data, where available. The tem- 
perature/luminosity cuts avoid complexities associated with 
variations in the fraction of baryonic matter expelled from 
the central regions of the clusters during their formation 
(Eke et al. 1998; Bialek, Evrard & Mohr 2001; it should 
be possible to relax these cuts in future work, given an im- 
proved calibration of the effect from simulations) . No quan- 
titative morphological classification procedure was used in 
the selection of the sample (although the inclusion of such 
a procedure using e.g. the power ratios method of Buote 
& Tsai 1995 would be straightforward in future work). Note 
that the selection function is not required for the determina- 
tion of cosmological parameters. We simply require accurate 
mass measurements. 

At moderate-to-high redshifts (2 > 0.3) the extension 
of the sample with respect to ASF02 was achieved primar- 
ily through two Large Programs of Chandra observations, 
lead by L. van Speybroeck and H. Ebeling, of clusters in 
the Massive Cluster Survey (MACS; Ebeling, Edge & Henry 
2001). From relatively short Chandra observations of 53 in- 
dividual MACS clusters, we identified 12 systems with a 
high degree of dynamical relaxation (details of the cluster 
morphologies are discussed by Ebeling et al. 2004, in prepa- 
ration). One of the clusters, MACS J 1423.8+2404, also has 
an additional deep Chandra observation which is used here. 
In addition to the MACS clusters, we have also included 
archival Chandra data for two other high redshift systems: 
MS1137.5+6625 (z = 0.782; Gioia & Luppino 1994) and 
C1J1226.9+3332 (2 = 0.892; Ebeling et al. 2001). The cen- 
tral X-ray emission in these clusters is less sharply peaked 
than most of the systems at lower redshifts (although the 
central cooling times in the clusters are still only 2 — 3 x 10 9 
yrs however). However, both clusters exhibit regular X-ray 
morphologies and are the most apparently relaxed, X-ray lu- 
minous clusters that we are aware of at such high redshifts. 
Additional support for the inclusion of MS1137. 5+6625 in 
our study comes from the agreement of the best fitting to- 
tal mass model determined from the Chandra X-ray data 
(presented here) and the independent weak lensing study 
of Clowe et al. (2000). For an SCDM cosmology, the Chan- 
dra data are well described by an NFW mass model with 
c = 3.5^j'g and r s = 185^7g 5 /i _1 kpc, implying r2oo = cr s = 
ttOtlf^h- 1 kpc. For the same cosmology, Clowe et al. (2000) 
find c = 4.2 and r2oo ~ 730h~ 1 kpc. The effective veloc- 
ity dispersion corresponding to the best-fit Chandra mass 
model, cr = 1100t|{$ km s 1 , is also consistent with the ob- 
served value of a = 8841^ kms -1 (Donahue et al. 1999). 
Although no detailed weak lensing study is available for 
C1J1226. 9+3332, Maughan et al. (2004) present a tempera- 
ture map from XMM-Newton observations which supports 
the identification of this system as a regular, dynamically re- 
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Table 1. Summary of the Chandra observations. The columns list the target name, observation date, detector used, observation mode 
and net exposure after all cleaning and screening processes were applied. The targets are listed in RA order. 





Date 


1 Intprt i~iv 

17CUCL L Ul 


Mode 




MACSJ0242.6-2132 


2002 Feb 07 


ACIS-I 


VFAINT 


10.2 


Abcll 383(1) 


2000 Nov 16 


ACIS-S 


VFAINT 


18.0 


Abcll 383(2) 


2000 Nov 16 


ACIS-I 


VFAINT 


17.2 


MACSJ0329.7-0212 


2002 Dec 12 


ACIS-I 


VFAINT 


16.8 


Abell 478 


2001 Jan 27 


ACIS-S 


FAINT 


39.9 


MACSJ0429.6-0253 


2002 Feb 07 


ACIS-I 


VFAINT 


19.1 


RXJ0439.0+0520 


2000 Aug 29 


ACIS-I 


FAINT 


7.6 


MACSJ0744.9+3927(1) 


2001 Nov 12 


ACIS-I 


VFAINT 


17.1 


MACSJ0744.9+3927(2) 


2003 Jan 04 


ACIS-I 


VFAINT 


14.5 


Abcll 611 


2001 Nov 03 


ACIS-S 


VFAINT 


34.5 


4C55 


2004 Jan 03 


ACIS-S 


VFAINT 


77.8 


MACSJ0947.2+7623 


2000 Oct 20 


ACIS-I 


VFAINT 


9.6 


Abell 963 


2000 Oct 11 


ACIS-S 


FAINT 


34.8 


MS1137.5+6625 


1999 Sep 30 


ACIS-I 


VFAINT 


103.8 


Abcll 1413 


2001 May 16 


ACIS-I 


VFAINT 


8.1 


C1J1226.9+3332 


2003 Jan 27 


ACIS-I 


VFAINT 


25.7 


MACSJ1311. 0-0311 


2002 Dec 15 


ACIS-I 


VFAINT 


12.0 


RXJ1347.5-1145(1) 


2000 Mar 03 


ACIS-S 


VFAINT 


8.6 


RXJ1347.5-1145(2) 


2000 Apr 29 


ACIS-S 


FAINT 


10.0 


RXJ1347.5-1145(3) 


2003 Sep 03 


ACIS-I 


VFAINT 


49.3 


Abell 1835(1) 


1999 Dec 11 


ACIS-S 


FAINT 


18.0 


Abell 1835(2) 


2000 Apr 29 


ACIS-S 


FAINT 


10.3 


3C295 


1999 Aug 30 


ACIS-S 


FAINT 


17.0 


MACSJ1423.8+2404 


2003 Aug 18 


ACIS-S 


VFAINT 


112.5 


Abell 2029 


2000 Apr 12 


ACIS-S 


FAINT 


19.2 


MACSJ1532.9+3021(1) 


2001 Aug 26 


ACIS-S 


VFAINT 


9.4 


MACSJ1532.9+3021(2) 


2001 Sep 06 


ACIS-I 


VFAINT 


9.2 


MACSJ1621.6+3810 


2002 Oct 18 


ACIS-I 


VFAINT 


7.9 


MACSJ1720.3+3536 


2002 Nov 03 


ACIS-I 


VFAINT 


16.6 


MACSJ1931. 8-2635 


2002 Oct 20 


ACIS-I 


VFAINT 


12.2 


MS2137.3-2353(1) 


1999 Nov 18 


ACIS-S 


VFAINT 


20.5 


MS2137.3-2353(2) 


2003 Nov 18 


ACIS-S 


VFAINT 


26.6 


MACSJ2229.8-2756 


2002 Nov 13 


ACIS-I 


VFAINT 


11.8 



cessed using the CI AO (version 3.0.2) software package, in- 
cluding the latest gain maps and calibration products. Bad 
pixels were removed and standard grade selections applied. 
Where possible, the extra information available in VFAINT 
mode was used to improve the rejection of cosmic ray events. 
For observations carried out with the ACIS-I detector, we 
have used the Chandra X-ray Center (CXC)/MIT charge 
transfer inefficiency correction. Time-dependent gain correc- 
tions were applied using A. Vikhlinin's apply_gain routine. 
The data were cleaned of periods of anomalously high back- 
ground using the same author's lc-dean script, using the 
recommended energy ranges and bin sizes for each detector. 
The net exposure times for the individual clusters are sum- 
marized in Table 1 and vary between 7.6 and 112.5 ks. The 
total good exposure time is 825.8 ks. 

The data have been analysed using the methods de- 
scribed by Allen et al. (2001a, 2002) and Schmidt et al. 
(2001). In brief, concentric annular spectra were extracted 
from the cleaned event lists, centred on the peaks of the X- 
ray emission from the clusters.* The spectra were analysed 

* For RXJ1347-1145, the data from the southeast quadrant of 



laxed cluster, well suited to the present work. The observed 
velocity dispersion of a — 997^205 kms -1 (Maughan et al. 
2004) is also consistent with the value of a > 1000 kms -1 in- 
ferred from the Chandra data. 

At low redshifts (z < 0.3) the extension of the sam- 
ple with respect to ASF02 was achieved via a search of the 
Chandra archive. Two clusters from the ASF02 sample have 
been dropped from this work: PKS0745-191 (z = 0.103) was 
dropped due to difficulties associated with extrapolating the 
observed profile beyond the restricted Chandra ACIS-S field 
of view. Abell 2390 (z = 0.230) was dropped due to dy- 
namical activity which is not localized and cannot easily 
be excluded from the analysis. The results for both clusters 
listed in ASF02 are, however, consistent with the analysis 
presented here. 



2.2 Chandra observations 

The Chandra observations were carried out using the Ad- 
vanced CCD Imaging Spectrometer (ACIS) between 1999 
August 30 and 2004 Jan 03. The standard level- 1 event lists 
produced by the Chandra pipeline processing were repro- 
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Table 2. The measured X-ray gas mass fractions at r2500 (68 
per cent confidence limits) for the reference SCDM and ACDM 
cosmologies. The results for Abell 2029 and 478 have been extrap- 
olated as described in the text. Redshifts for the MACS clusters 
are from Ebcling et al. (2004, in preparation). A table contain- 
ing the redshift and / gas data is available from the authors on 
request. 





SCDM 
f ft" 18 


ACDM 
f ft- 1 - 8 




0.175 


± 


0.023 


0.130 


± 


0.018 


/A_LH^11 ooo 


0.169 


± 


0.011 


0.122 


± 


0.009 


MACSTfn99 7-0919 


0.155 


± 


0.019 


0.119 


± 


0.018 


A hell 478 


0.184 


± 


0.011 


0.120 


± 


0.008 


MACS1042Q 6-0253 


0.177 


± 


0.017 


0.141 


± 


0.015 


RX 10439 0+0520 


0.137 


± 


0.018 


0.098 


± 


0.013 


MACS 10744 9+3927 


0.155 


± 


0.018 


0.141 


± 


0.019 


Ahell fi1 1 


0.149 


± 


0.017 


0.111 


± 


0.016 


4C55 


0.163 


± 


0.009 


0.115 


± 


0.005 


MACS 10947 2+7623 


0.173 


± 


0.019 


0.130 


± 


0.016 


AK__,11 QfiQ 

i-vucii yoo 


0.180 


± 


0.015 


0.128 


± 


0.012 


MS1137.5+6625 


0.100 


± 


0.016 


0.094 


± 


0.014 


Abell 1413 


0.167 


± 


0.019 


0.114 


± 


0.013 


C1J1226. 9+3332 


0.114 


± 


0.021 


0.102 


± 


0.027 


MACSJ1311. 0-0311 


0.094 


± 


0.025 


0.072 


± 


0.022 


RXJ1347.5-1145 


0.137 


± 


0.009 


0.108 


± 


0.009 


Abell 1835 


0.164 


± 


0.012 


0.112 


± 


0.012 


3C295 


0.129 


± 


0.019 


0.106 


± 


0.018 


MACSJ1423.8+2404 


0.135 


± 


0.011 


0.113 


± 


0.008 


Abell 2029 


0.189 


± 


0.011 


0.121 


± 


0.008 


MACSJ1532.9+3021 


0.159 


± 


0.017 


0.114 


± 


0.017 


MACSJ1621. 6+3810 


0.156 


± 


0.034 


0.131 


± 


0.029 


MACSJ1720.3+3536 


0.159 


± 


0.024 


0.123 


± 


0.020 


MACSJ1931. 8-2635 


0.189 


± 


0.025 


0.145 


± 


0.022 


MS2137.3-2353 


0.169 


± 


0.010 


0.124 


± 


0.009 


MACSJ2229.8-2756 


0.177 


± 


0.018 


0.139 


± 


0.017 



using XSPEC (version 11.3: Arnaud 1996), the MEKAL 
plasma emission code (Kaastra & Mewe 1993; incorporat- 
ing the Fe-L calculations of Liedhal, Osterheld & Gold- 
stein 1995) and the photoelectric absorption models of 
Balucinska-Church & McCammon (1992). The ACISABS 
model was used to account for time-dependent contamina- 
tion along the instrument light path. We have allowed for the 
small amount of extra contamination present in ACIS-I ob- 
servations by increasing the optical depth of the ACISABS 
model contaminant by At = 0.14 at 0.67 keV (A. Vikhlinin, 
private communication). Only data in the 0.8 — 7.0 keV en- 
ergy range were used for our analysis (the exceptions being 
the earliest observations of 3C295, Abell 1835 and Abell 
2029 where a wider 0.6-7.0 keV band was used). The spec- 
tra for all annuli were modelled simultaneously, in order to 
determine the deprojected X-ray gas temperature profiles 
under the assumption of spherical symmetry. 

For the nearer clusters (z < 0.3), background spectra 
were extracted from the blank-field data sets produced by M. 
Markevitch and available from the CXC. For the more dis- 
tant systems (and the first observation of Abell 1835, which 



the cluster were excluded due to ongoing merger activity in that 
region (see Allen et al. 2002). 



has an unusually high, but relatively constant, background 
level) background spectra were extracted from appropriate, 
source free regions of the target data sets. (We have con- 
firmed that similar results are obtained using the blank- 
field background data sets throughout.) Separate photon- 
weighted response matrices and effective area files were con- 
structed for each region studied, using the calibration files 
appropriate for the period of observation. For the ACIS-I 
analysis, we have decreased the quantum efficiency at ener- 
gies below 2 keV by 7 per cent from the nominal value, and 
then at all energies by a further 8 per cent. This improves the 
cross-calibration between ACIS-S and ACIS-I observations 
of clusters in our sample and is consistent with the recom- 
mendations of the Chandra calibration team (A. Vikhlinin, 
private communication) . 



2.3 Chandra / gaa measurements 

For the mass modelling, azimuthally-averaged surface 
brightness profiles were constructed from background sub- 
tracted, flat-fielded images with a 0.984 x 0.984 arcsec 2 pixel 
scale (2x2 raw detector pixels) . When combined with the de- 
projected spectral temperature profiles, the surface bright- 
ness profiles can be used to determine the X-ray gas mass 
profiles (to high precision) and total mass profiles in the 
clusters. ^ For this analysis, we have used an enhanced ver- 
sion of the image deprojection code described by White, 
Jones & Forman (1997). 

We have parameterized the cluster total mass profiles 
(luminous plus dark matter) using a Navarro, Frenk & White 
(1997; hereafter NFW) model with 



p(r) = 



pc(z)S c 



(r/r s ) (1 + r/r s y 



(1) 



where p(r) is the mass density, p c (z) = 3H(z) 2 /&ivG is the 
critical density for closure at redshift z, r s is the scale radius, 
c is the concentration parameter (with c = T2oo/r B ) and 
S c = 200c 3 /3 [ln(l + c) - c/(l + c)]. The normalizations of 
the mass profiles may also be usefully expressed in terms of 
an effective velocity dispersion, a — \/50r B cH(z) (with r B 
in units of Mpc and H(z) in kms" 1 Mpc -1 ). Mass models 
were examined over regular 100 x 100 grids in the (r s ,a) 
plane. 

In determining the results on the X-ray gas mass frac- 
tion, /g as , we have adopted a canonical radius of r2soo, within 



T The observed surface brightness profile and parameterized mass 
model arc together used to predict the temperature profile of the 
X-ray gas. We use the median temperature profile determined 
from 100 Monte-Carlo simulations. The outermost pressure is 
fixed using an iterative technique which ensures a smooth pres- 
sure gradient in these regions. The predicted temperature profile 
is rebinned to the same binning as the spectral results and the \ 2 
difference between the observed and predicted, deprojected tem- 
perature profiles is calculated. The parameters for the mass model 
are stepped through a regular grid of values in the r B -a plane 
(see text) to determine the best-fit values and 68 per cent confi- 
dence limits. The gas mass profile is determined to high precision 
at each grid point directly from the observed surface brightness 
profile and model temperature profile. Spherical symmetry and 
hydrostatic equilibrium arc assumed throughout. 
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Figure 2. The apparent variation of the X-ray gas mass fraction (with root- mean-square la errors) as a function of redshift for the 
reference (a: left panel) SCDM and (b: right panel) ACDM cosmologies. The grey curve in (a) shows the predicted / gas (z) behaviour 
for the best-fitting model cosmology with f2 m = 0.25 and = 0.96 (see Section 3.2). Clusters at higher redshifts appear to have lower 
gas mass fractions because the SCDM cosmology underestimates the relative distances to these systems. The curve in (b) shows the 
best-fitting constant value. 



o 




0.5 1 

r / r 2500 



Figure 1. The observed X-ray gas mass fraction profiles for the 26 
clusters with the radial axes scaled in units of r2soo ■ The reference 
ACDM cosmology is assumed. Note that / gas (r) is an integrated 
quantity and the error bars on neighbouring points in a profile 
are correlated. 



which the mean mass density is 2500 times the critical den- 
sity of the Universe at the redshifts of the clusters. The 
f 2500 values are determined directly from the Chandra data, 
with confidence limits calculated from the \ 2 grids and are 
(in general) well-matched to the outermost radii at which 
reliable temperature measurements can be made from the 
Chandra data. 



Fig. 1 shows the observed / gas (r) profiles for the 26 clus- 
ters in our sample, determined from the Chandra data using 
the reference ACDM cosmology. Although some variation is 
present from cluster to cluster, particularly at small radii, 
the profiles tend towards a universal value at r2soo- (We 
note that the cluster with the most discrepant f sas (r) profile 
and highest / gas value at small radii is MACSJ1532. 9+3021, 
which exhibits an unusually high ellipticity in its innermost 
regions. 4C55 also exhibits a high ellipticity and isophote 
shifts at small radii and has the second highest / gas value at 
0.1f"25oo. Both clusters appear relaxed at larger radii, how- 
ever, and their / gas profiles have recovered to the 'universal' 
form by T2500-) Table 2 lists the results on the X-ray gas 
mass fractions at r25oo for the reference SCDM and ACDM 
cosmologies. Taking the weighted mean of the ACDM results 
we obtain / gas = 0.1173 ± 0.0022 /if 15 . 

In calculating the total baryonic mass in the clusters, we 
assume that the optically luminous baryonic mass in galax- 
ies scales as 0.19/i ' 5 times the X-ray gas mass. This result 
is based on detailed studies of nearby and intermediate red- 
shift clusters (White et al. 1993, Fukugita, Hogan & Peebles 
1998; see also Voevodkin & Vikhlinin 2004) and corresponds 
to ~ 16 per cent of the X-ray gas mass. Uncertainties in this 
correction have a negligible impact on the overall error bud- 
get. Other sources of baryonic matter in the clusters are 
expected to make very small contributions to the total mass 
and are ignored. 

We note that the Chandra data for Abell 478 and 2029 
do not extend quite to r^soo- For these clusters, we measure 
/ gaa directly at r = 0.75 r2soo for the SCDM cosmology or 
r = 0.70 r25oo for the ACDM cosmology and extrapolate the 
results to r2soo using the median / gas (7*) profile determined 
from Fig. 1. This extrapolation results in corrections to the 
directly measured / gas values of ~ 5 per cent. To be conser- 
vative, we have included a 5 per cent systematic uncertainty 



© 0000 RAS, MNRAS 000, 000-000 



6 S. W. Allen et al. 



in the tabulated / gas measurements for Abell 478 and 2029 
to allow for uncertainties in this extrapolation. 

Fig. 2 shows the / gas values as a function of redshift 
for the reference SCDM and ACDM cosmologies. Whereas 
the results for the ACDM cosmology are consistent with a 
constant / gas value (\ 2 = 22.7 for 25 degrees of freedom) 
the results for the reference SCDM cosmology indicate an 
apparent drop in / gas as the redshift increases. The \ 2 = 
61.8 obtained from a fit to the SCDM data with a constant 
model indicates that the SCDM cosmology is inconsistent 
with the expectation that / gas (z) should be constant. 



3 COSMOLOGICAL ANALYSIS 
3.1 Dark energy models 

We have considered two separate dark energy models in our 
analysis: standard ACDM and the extended XCDM param- 
eterization. Our definitions of the relevant quantities closely 
follow Peebles & Ratra (2003) and we refer the reader to 
that work for a discussion of the underlying assumptions. 

The Friedmann equation, which relates the first 
time derivative of the scale factor of the Universe, a, 
to the total density can be conveniently expressed as 
(a/a) 2 = H(z) 2 = HqE(z) 2 , where 



E(z) = ^JVL m (l + z) 3 + fi x /(z) + fi K (l + z) 



(2) 



2.4 CMB analysis 

Our analysis of CMB observations uses the WMAP temper- 
ature (TT) data for multipolcs I < 900 (Hinshaw et al. 2003) 
and temperature-polarization (TE) data for / < 450 (Kogut 
et al. 2003). To extend the analysis to higher multipolcs 
(smaller scales), we also include data from the Cosmic Back- 
ground Imager (CBI; Pearson et al. 2003) and Arcminute 
Cosmology Bolometer Array Receiver (ACBAR; Kuo et al. 
2003) for / > 800. The comparison of model angular power 
spectra with the WMAP data employs the likelihood calcu- 
lation routines released by the WMAP team (Verde et al. 
2003). 

Our analysis of the CMB data uses the CosmoMC 
code'f . This in turn uses CAMB (Lewis, Challinor & Lasenby 
2000), which is based on CMBFAST (Seljak & Zaldar- 
riaga 1996), to generate the CMB and matter power spec- 
trum transfer functions, and a Metropolis-Hastings Markov 
Chain Monte Carlo (MCMC) algorithm to explore parame- 
ter space. We used the covariance matrix of the parameters 
calculated from an initial set of test runs to improve sam- 
pling efficiency (see Lewis & Bridle 2002 for more details). 

We have fitted the data using an extended XCDM cos- 
mological model with eight free parameters: the physical 
dark matter and baryon densities in units of the critical 
density, the curvature fix, the Hubble constant, the dark en- 
ergy equation of state parameter, the recombination redshift 
(at which the reionization fraction is a half, assuming in- 
stantaneous reionization) , the amplitude of the scalar power 
spectrum and the scalar spectral index. We also examined 
a standard ACDM model in which we fixed i»x = — 1. In 
all cases, we have assumed an absence of tensor components 
and included uniform priors 30 < Hokms" 1 Mpc _1 < 100 
and —4 < wx < 1. (Tests in which tensor components were 
included with ACDM models lead to similar results on dark 
energy, but took much longer to compute.) 

The analysis was carried out on the Cambridge X-ray 
group Linux cluster. For each model we accumulated a total 
of at least 10 6 correlated samples in 10 separate chains. We 
satisfied ourselves that the chains had converged by ensur- 
ing that consistent final results were obtained from numer- 
ous small subsets of the chains. In all cases, we allowed a 
conservative burn-in period of 10 4 samples for each chain. 



i http://cosmologist.info/cosmomc/ 



Here fix is the dark energy density and f(z) its redshift 
dependence. (We have ignored the density contribution from 
radiation and relativistic matter.) For ACDM cosmologies, 
the dark energy density is constant and f(z) = 1. Within the 
extended XCDM dark energy parametrization, the pressure 
is related to the density as px = wxpx so that for constant 
wx, the dark energy density scales as px oc a -3 ^""*) and 
f(z) = (l + z) 3(1+u,x) . We note that for wx < -1/3 the dark 
energy makes a positive contribution to the acceleration of 
the expansion of the universe. For wx < — 1 the dark energy 
density is increasing with time. 

Our analysis of Chandra / gas data requires the angular 
diameter distances to the clusters which are defined as 



d,A = 



H (l + z) 



sinh 



(3) 



3.2 Analysis of the / gas data 

The differences between the shapes of the / gas (z) curves in 
Figs. 2(a) and (b) reflect the dependence of the measured 
/gas (2) values on the angular diameter distances to the clus- 
ters: / gas oc (1 + 2) 2 d A ' 5 . Under the assumption (Section 1) 
that / gas should, in reality, be constant with redshift, sim- 
ple inspection of Fig. 2 clearly favours the ACDM over the 
SCDM cosmology. 

To determine constraints on the relevant cosmological 
parameters, we have fitted the f sas (z) data in Fig. 2(a) with 
a model that accounts for the expected apparent variation 
in / gas (z) as the underlying cosmology is varied. (We choose 
to work with the SCDM data as our reference cosmology, 
although similar results can be derived using the ACDM 
data set.) Note that the f sas (r) profiles exhibit only small 
variations around T2500 so changes in r25oo as the cosmology 
is varied can be ignored. The model function fitted to the 
data is 



f!^ M (z) = 



b fi b 



'l + 0.19V7i) fi„ 



.SCDM 
"A 



Jmod 
"A 



(4) 



where d^ od (z) and d A CDM (2) are the angular diameter dis- 
tances to the clusters in the current model and reference 
SCDM (h = 0.5) cosmologies. Note that although varia- 
tions in the dark energy density affect only the shape of the 
/gas(z) curve, the normalization depends on fi m , fib, h and 
b, where b is a bias factor motivated by gasdynamical sim- 
ulations which suggest that the baryon fraction in clusters 
is slightly lower than for the universe as a whole (e.g. Eke, 
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Figure 3. The X-ray bias factor, b (the enclosed baryon fraction 
relative to the universal value) as a function of radius in units of 
the virial radius r v ; r , from the simulations of Eke et al. (1998). 
The simulated clusters have similar masses to the systems studied 
here. The results for the three most dynamically relaxed clusters 
in the simulations are shown as darker curves. The solid circles 
mark the median profile determined from the Chandra observa- 
tions, scaled from Fig. 1 assuming Q m = 0.25, H b = 0.0413 and 
r 2500 = 0.25 r v ; r . Beyond a radius r > 0.2r v j r , the simulated clus- 
ters exhibit consistent, relatively flat b profiles. At r = 0.25r v j r , 
a radius comparable to the measurement radius of the Chandra 
observations, the simulations give b = 0.824 ± 0.033. 



Navarro & Frenk 1998; Bialek et al. 2001). We use the re- 
sults of Eke et al. (1998) from simulations of 10 clusters of 
similar masses to the observed systems to constrain 6. Ex- 
cluding the data for the most dynamically active cluster in 
that study (recall that the / gas data are drawn from Chandra 
observations of dynamically relaxed systems) , the simulated 
clusters show consistent, relatively flat baryonic mass frac- 
tion profiles for radii r > 0.2r v i r (Fig 3). At r = 0.25r v i r , a 
radius comparable to the measurement radius for the Chan- 
dra observations, the simulations of Eke et al. (1998) give 
6 = 0.824 ± 0.033. 

We note the excellent agreement between the median, 
scaled f ga s(r) profile determined from the Chandra data 
(shown as dark circles in Fig 3) and the simulated profiles 
for the three most relaxed clusters in the study of Eke et al. 
(1998; the darker curves in Fig 3). This agreement supports 
the use of the simulations in estimating the bias factor at 
f2500- Note also that the simulations of Eke et al. (1998) in- 
dicate negligible evolution of the bias parameter (measured 
within r ~ 0.5r v i r ) over the redshift range considered here. 

For our analysis of the Chandra / gaa data alone, we 
employ simple Gaussian priors on fib ft 2 and h. Two sepa- 
rate sets of priors were used: 'standard' priors with fib ft 2 = 
0.0214 ± 0.0020 (Kirkman et al. 2003) and h = 0.72 ± 0.08 
(Freedman et al. 2001), and 'weak' priors in which we triple 
the nominal uncertainties to give Q b h 2 = 0.0214 ± 0.0060 



and h = 0.72 ± 0.24. (No priors on fib/i 2 and h are assumed 
in the combined analysis of / gas and CMB data; Section 3.3.) 

We assume a Gaussian prior on b. The rms fractional 
deviation in b from the simulations of Eke et al. (1998; ~ 4 
per cent) was added in quadrature to a nominal 10 per cent 
systematic uncertainty associated with the overall normal- 
ization of the /gas («) curve. This allows for residual uncer- 
tainties associated with the simulations and/or the calibra- 
tion of the Chandra instruments. § Thus, for the analysis of 
the f s 

for any particular model is 



data alone using the standard priors, the x 2 value 



+ 



^ [/IS DM (*i) 



fi b /i 2 - 0.0214 \ 2 Z/i-0.72 



0.0020 



+ 



0.08 



-)' + ( 



b - 0.824 
0.089 



(5) 



Here / ga s, i and 07 aa ; are the observed values and symmet- 
ric rms errors for the SCDM cosmology from Table 2. Note 
that we have not accounted for the intrinsic cluster-cluster 
scatter in b explicitly. However, the simulations suggest that 
this scatter is small when compared with the statistical un- 
certainties in the / gas measurements. 



3.3 Combination of / gas and CMB constraints 

For the combined Chandra+CMB analysis, we importance 
sample the MCMC results from the CMB analysis, folding 
in the / gas constraints (Allen, Schmidt & Bridle 2003). Each 
of the MCMC samples from the CMB analysis provides a 
value for fi m , fix, wx, Ho and fib ft 2 . Using these values, 
we fit the fg as (z) data with the model described by Equa- 
tion 4, including the same Gaussian prior on the bias factor, 
including the allowance for systematic uncertainties in the 
normalization of the curve. This provides a x 2 value for each 
of the MCMC samples. The weight of the MCMC sample is 
then multiplied by e x ' . 



4 COSMOLOGICAL CONSTRAINTS 

4.1 Results for the ACDM cosmology 

For the ACDM cosmology, we have examined a grid of cos- 
mological models covering the plane 0.0 < fi m < TO and 
0.0 < fi A < 2.0. The joint 68.3, 95.4 and 99.7 per cent con- 
fidence contours (corresponding to Ax 2 values of 2.30, 6.17 
and 11.8, respectively) obtained from the Chandra / gas data, 
including standard priors on fib ft 2 (Kirkman et al. 2003) and 
h (Freedman et al. 2001) are shown in Fig. 4. The best-fit 
parameters and marginalized 68 per cent confidence limits 



obtained using the standard priors are fi n 



0.245 



+ 0.040 



3 The agreement between the independent mass measurements 
from X-ray and gravitational lensing studies for several of the 
target clusters argues that the systematic uncertainties are un- 
likely to significantly exceed 10 per cent. 
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Clusters 

(+BBNS+HST) 




Figure 4. The 68.3, 95.4 and 99.7 per cent (1, 2 and 3tr) confi- 
dence constraints in the f2 m , Q\ plane obtained from the analysis 
of the cluster / gaa data using standard priors on Q^h 2 (Kirkman 
et al. 2003) and h (Freedman et al. 2001). Also shown are the in- 
dependent results obtained from CMB data using a weak uniform 
prior on h (0.3 < h < 1), and Type la supernovae data (Tonry 
et al. 2003). A ACDM cosmology is assumed with the curvature, 
Ok; included as a free parameter in the analysis. 



Figure 6. The 68.3, 95.4 and 99.7 per cent confidence con- 
straints in the Q m ,Q\ plane obtained from the analysis of the 
combined / gas +CMB data set using the ACDM model. We find 
marginalized 68 per cent confidence limits of f2 m = 0.28 
Q A = 0.73l u 'o3, with U K = -0.01 ± 0.02. 



+ 05 and 
0.04 ana 



X 




Figure 5. The marginalized constraints on Qa obtained from 
the Chandra / gas data using the standard (dark curve) and weak 
(light curve) priors on Q^h 2 and h. The dotted lines mark the 
formal 1, 2 and 3<r limits (see text for details). The curvature is 
a free parameter in the analysis. 



X 



to 




Figure 7. The marginalized constraints on !!k (= 1 — O m — S1a) 
obtained from the Chandra / gas data using the standard priors 
on Q b /i 2 and h. The dotted lines mark the formal 1, 2 and 3<r 
limits (see text for details). 
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and A = 0.96t£g, with x 2 = 24.5 for 24 degrees of free- 
dom. The x 2 value indicates that the model provides an 
acceptable description of the data. 

Fig. 5 shows the marginalized constraints on Oa ob- 
tained using both the standard and weak priors on Ob ft 2 
and h. We see that even using the weak priors (Ob/i 2 = 
0.0214 ± 0.0060, h = 0.72 ± 0.24), the / gas data still provide 
a clear detection of the effects of Oa at > 3cr significance 



0.94 



+0.21^ 



A Monte Carlo analysis of the data indi- 



cates that Oa < is ruled out at > 99.9 per cent confidence. 

Fig. 4 also shows a comparison with independent con- 
straints obtained from the CMB data using only a weak 
uniform prior on h (0.3 < h < 1.0), and from Type la super- 
novae studies (Tonry et al. 2003). The agreement between 
the / gas and CMB constraints in particular is reassuring and 
motivates the combined analysis of these data sets, discussed 
below. 

Fig. 6 shows the constraints on O m and Oa obtained 
from the analysis of the combined / gas +CMB data set. No 
priors, other than the constraint b — 0.824 ± 0.089 are as- 
sumed. We see that the / gas +CMB data set provides a re- 
markably tight constraint in the O m ,OA plane, with best 
fit values and marginalized 68 per cent confidence limits of 
O m = 0.28t°-°ij and A = 0.73±°;°*. 

Fig. 7 shows the marginalized constraints on Ok (= 
1 — O m — Qa) from the / gas data using the standard priors 
on fl b h 2 and h. The best fit result of K = -0.2 ± 0.2 
(68 per cent confidence limits) is consistent with the much 
tighter constraint of Ok = —0.01 ± 0.02 obtained from the 
combined /gas+CMB data set. (The CMB data alone give 



-0.03 i 



using only a wide uniform prior on the 



Hubble constant). 

Finally, Fig. 8 shows an alternative way to visualize the 
effects of dark energy on the distances to the clusters. The 
model function shown in the figure is dA° d ( z )//gas d ( z = 



0) 2/3 and the data points d s A CDM {z)/ f°* B {z) 2/3 . The use of 
this function removes all dependence on the reference cos- 
mology in the figure. The dark curve in Fig. 8 shows the 
results for the best-fitting ACDM cosmology. Also shown 
are the curves obtained keeping the parameters fixed at 
their best-fit values, but setting Oa = (grey curve; this 
shows the effects of the dark energy component), and set- 
ting Oa = 0, Q m = 1 (dotted curve). The model includ- 
ing the dark energy component provides the best descrip- 
tion of the data over the full redshift range. The Oa = 0, 
O m = 1 model provides a poor fit, and would require 
Ho ~ 20kms _1 Mpc -1 to approximately match the ob- 
served normalization (given the standard set of priors). 



4.2 Extended XCDM models 

Fig. 9 shows the 68.3, 95.4 and 99.7 per cent confidence con- 
straints in the O m , Ox plane for the extended XCDM models 
from the analysis of the combined / gas +CMB data set. We 
obtain best fitting values and marginalized 68 per cent con- 
fidence limits of Ox = 0.75 ± 0.04 and O m = 0.26^ The 
constraints on the mean matter and dark energy densities for 
the extended XCDM models are similar to those obtained 
for the ACDM cosmology. The curvature is measured to be 
K = -0.02 ±0.02. 

Fig. 10(a) shows the constraints in the O m , wx plane ob- 




Figure 8. (Upper panel) The function d^ od (z)/[/™° d (^ = 0)] 2 / 3 
for the best-fitting ACDM cosmology (dark curve). Also shown 
are the curves obtained keeping all parameters fixed at their best- 
fit values but setting Qa = (grey curve), and setting Qa = 0, 
Qm = 1 (dotted curve). (Lower panel) The ratios with respect to 
the best-fitting ACDM model. 




Figure 9. The 68.3, 95.4 and 99.7 per cent confidence constraints 
in the Q m , Qx plane obtained from the analysis of the combined 
/gas+CMB data set. We obtained marginalized 68 per cent confi- 



dence limits of Q n 



: 0.26 



+ 0.06 



and Q x = 0.75±0.04. The results 



are similar to those obtained for the ACDM models in Fie 
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Figure 10. The 68.3, 95.4 and 99.7 per cent confidence constraints in the £7 m ,ui x plane from the analysis of (a: left panel) the combined 
/gas+CMB data, with !2j; free. Also shown, for comparison purposes, are the results obtained from the analysis of Type la supernovae 
data by Tonry et al. (2003) assuming a flat geometry, (b: right panel) The results obtained from the Chandra / gas data alone, assuming 
standard priors on Cl^h 2 and h and a flat geometry, together with the supernovae constraints. 




Figure 11. (a: left panel) The marginalized constraints on the equation of state parameter, t«x, obtained from the / gas +CMB data 
with Ck f ree (dark solid curve), / gas +CMB data assuming a flat geometry (dotted curve) and / gas data alone assuming a flat geometry 
and standard priors on Q^h 2 and h (grey curve), (b) As for (a) but imposing the prior constraint wx > — 1. 
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tained from from the same data. Also shown, for comparison 
purposes, are the results from Type la supernovae studies 
(Tonry et al. 2003). Fig. 10(b) shows the results obtained 
from the / gas data alone, assuming a flat geometry and the 
standard priors on fit,/) 2 and h. The results are in excellent 
agreement with those obtained from the / gas +CMB data 
set. We again note the ability of the / ga a data, used in com- 
bination with the CMB data or standard priors, to break 
important degeneracies between parameters. 

Fig. 11(a) shows the marginalized constraints on i»x for 
the extended XCDM models. For the / gas +CMB data with 
fix free we find wx = — 1.26to'.24- Under the assumption of 
a flat geometry, the same / gas +CMB data data give wx = 
— 1.221q'22- For a flat geometry, the / gaa data and standard 
priors on fib/i 2 and h give wx = — 1. 20^28- Note that for a 
flat geometry, the supernovae data alone give w x = —2. 

21?:! 

(68 per cent confidence limits). 

Fig. 11(b) shows the marginalized constraints on wx 
obtained from XCDM models when we apply the prior con- 
straint wx > — 1. Under this assumption, the / gas -|-CMB 
data with fix free give a 95.4 per cent confidence constraint 
of wx < —0.69. If we assume flatness, the same data require 
wx < —0.75. For a flat geometry and standard priors on 
fib/i 2 and h, the / gas data alone give tux < —0.69. These con- 
straints are similar to those obtained by Tonry et al. (2003; 
wx < —0.73) from supernovae data using a prior on fi m 
from the 2dF Galaxy Redshift Survey (fi m fr = 0.20 ± 0.03; 
Percival et al. 2001) and assuming a flat geometry. Our re- 
sults on wx are also consistent with (and comparable to) 
those reported by the WMAP team (Spergel et al. 2003). 

Finally, we note that our results for the extended 
XCDM cosmology imply that the mean matter and 
dark energy densities become equal at a redshift z = 
(On/fix) 1731 "* - 1 = 0.30 ± 0.09, and that the Universe 
moves from a decelerating to an accelerating phase at z = 
[-(1 + 3m x )!ix/fi m r 1/3, " x - 1 = 0.70 ± 0.11 (68 per cent 
confidence limits). 



5 DISCUSSION 

Our results provide the first clear confirmation of type 
la supernovae results in terms of detecting the effects of 
dark energy on distance measurements to a separate, well- 
defined source population. Our results cover the redshift 
range where the expansion of the Universe moves from a 
decelerating to an accelerating phase. The significance of 
our detection of dark energy is > 3a (> 99.9 per cent sig- 
nificance from Monte Carlo simulations) for the standard 
ACDM model with fix free, using only weak priors on Qbh 2 
and h. This accuracy is comparable to that obtained from 
current supernovae work (Tonry et al. 2003; see also Riess 
et al. 2004). 

It is interesting to note that our preferred value for tux 
in the extended XCDM models is slightly less than -1, which 
allows the possibility that the dark energy density is increas- 
ing with time. Such a scenario is also mildly favoured by 
recent Type la supernovae studies (Tonry et al. 2003; Riess 
et al. 2004). We stress, however, that the Chandra results 
remain consistent with ACDM (wx = — 1). 

A major benefit of our technique is that the applica- 
tion of standard priors on fib ft 2 and Ho, or the combina- 



tion with CMB data, also leads to tight constraints on fi m , 
thereby allowing important degeneracies between parame- 
ters to be broken. For a ACDM cosmology (fix free), we 
find fi m = 0.245ly;o37 using standard priors on fib/i 2 and h, 
or fi m = 0.28lg'o4 when the / gas and CMB data are com- 
bined. These constraints are comparable to those obtained 
from the combination of current CMB data with a variety 
of other data sets and priors (e.g. Spergel et al. 2003). We 
note that the lower value of fi m obtained in this work with 
respect to ASF02 is primarily due to the inclusion of the 
bias factor b in the present study, together with changes in 
the prior on Q b h 2 . The (slightly) larger error bars on fi m 
are due to the inclusion of the 10 per cent systematic uncer- 
tainty in the normalization of the / gas curve (via b), which 
is motivated by residual uncertainties in the calibration of 
the Chandra detectors. It may be that this 10 per cent al- 
lowance overestimates the systematics errors. If it were not 
included, the constraint on fi m for the ACDM cosmology 
from the / gas data using standard priors on Qbh 2 and h 
would become fi m = 0.246^0 029- Recall that the constraint 
on fi m arises primarily from the normalization of the / gas (z) 
curve and so is affected by the 10 per cent systematic un- 
certainty, whereas the constraint on fix is determined by 
the shape of the curve and is so largely independent of the 
uncertainties in Qbh 2 , h and b. 

The evidence for dark energy from the Chandra data 
is robust against uncertainties in the bias factor, b. The re- 
sults depend primarily upon the shape of the / gas (z) curve 
and so doubling the overall uncertainty in b to 20 per cent 
has little effect. Only redshift evolution in b can change the 
results on dark energy. However, to remove the evidence for 
dark energy (i.e. measure OX — 0) we would require 6 to 
decrease with increasing redshift by > 30 per cent over the 
interval < z < 1. This change in 6 is much larger than 
is allowed by simulations; the study of Eke et al. (1998) in- 
dicates negligible evolution over the redshift range studied 
here. For illustration purposes only, we have examined the 
effects of including (substantial) evolution in 6 such that 
b(z) = (1 — O.lz)b(O). This leads to only a small change in 
the results: fiA = 0.72^27 for the ACDM cosmology using 
the Chandra data and weak priors. (The detection of dark 
energy remains significant at the ~ 2.5 sigma level.) Includ- 
ing such evolution in 6 also shifts the best-fit value for w 
closer to -1 : w = — 0.98to' 2 4 for the same data using the 
XCDM model. 

An important aspect of the present work is that the 
clusters studied are regular, apparently dynamically relaxed 
systems. This results in a significant reduction of the scatter 
in the / gas measurements with respect to studies that do not 
include such a selection criterion (e.g. Ettori et al. 2003). 
Note also that our analysis does not impose a parametric 
form for the X-ray gas distribution, uses a realistic param- 
eterization for the total matter distribution, and makes full 
use of information on the temperature profiles in the clus- 
ters. Independent confirmation of the total masses within 
r'2500 is available from weak lensing studies in a number of 
cases, which lends support to the reliability of the / gas mea- 
surements (see discussion in ASF02; a program to expand 
the weak lensing measurements to the entire sample stud- 
ied here is underway.) Finally, we note that the effects of 
departures from spherical symmetry on the / gas results are 
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expected to be small (< a few per cent; Buote & Canizares 
1996, Piffaretti, Jetzer & Schindler 2003). 

An ASCII table containing the redshift and / gas (z) data 
is available from the authors on request. 
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